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The KTeV collaboration at Fermi National Accelerator Laboratory has recently completed 
searches for and measurements of several decay modes of the neutral kaon and pion. These 
include new searches for lepton flavor violating decays (which have not been seen), and a new 
study of the parity properties of the decay 7r° — > e + e~e + e~ . 



1 The KTeV Detector 



Fermilab's KTeV detector (Fig. [T]) was constructed for Experiments 799 and 832. The two 
experiments were designed to concentrate on different aspects of neutral kaon physics: E799 
on rare decays of the Kl and E832 on measurement of Re(e'/e). A primary proton beam with 
energy 800 GeV struck a BeO target at a targeting angle of 4.8 mrad, and collimation and 
sweeping magnets produced two parallel neutral hadron beams. The beams entered a 60 m long 
vacuum decay region, which ended at a Mylar-Kevlar vacuum window. Decay products were 
tracked with a series of drift chambers surrounding a dipole analysis magnet. Downstream of the 
drift chambers were a series of transition radiation detectors (TRD) (in E799 only) and a pure 
Csl electromagnetic calorimeter, an acive hadron beam absorber, and a set of muon detectors 
behind steel shielding. Photon veto detectors surrounded the fiducial volume in the transverse 
directions. The detector is described in more detail in Ref.H 



2 The decay ir — ► e + e e + e and the parity of the ir 

The neutral pion's parity has his toric ally been studied in two ways: indirectly via the cross- 
section of ir~ capture on deuteronJ^^, or directly via the double Dalitz decay ir° — > e + e~e + e~®. 
While both sets of results are consistent with the negative parity, the direct measurement has 
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Figure 1: The KTeV spectrometer as configured for E799. 
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Figure 2: Lowest order Feynman diagram for ir — > e e e e . The direct contribution is shown; a second 

diagram exists with e\ and exchanged. 

only 3.6 a significance. KTeV has now reported results^ that conclusively confirm the neg- 
ative 7T° parity as well as the first-ever searches for parity and CPT violaton, and the first 
measurements of the electromagnetic form factor, in this mode. 

The 7T° — > e + e~e + e~ decay proceeds through a two-photon intermediate state (Fig. [2]). The 
most general interaction Lagrangian for the ir° — > 7*7* transition can be written^H 

C oc C^F^F^ (1) 

where F^ v and F pa are the photon fields, $ is the pion field, and the coupling has the form 

C^pa oc /(xi, x 2 )[cos C^fiupa + sin (e lS {g w g va - g^gup)]- (2) 

The first term in C^pa is the expected pseudoscalar coupling and the second term introduces a 
scalar coupling with a mixing angle C an d a phase difference 5. Nuclear parity violation would 
introduce a nonzero while CPT violation would cause the phase 8 to be nonzero. We assume 
the standard parity-conserving form for the 7* — > e + e~ conversion. 

The form factor f(xi,X2) is expressed in terms of the momentum transfer of each of the 
virtual photons, or equivalently the invariant masses of the two Dalitz pairs: x\ = (m e + e - /M n o) 2 ; 

X2 = (m + -/M 7T o) 2 . In calculating the phase space variables for an individual event, there is 
an intrinsic ambiguity in assigning each electron to a positron to form a Dalitz pair. KTeV's 



analysis uses a matrix element model that includes the exchange diagrams and therefore avoids 
the need to enforce a pairing choice. The form factor is parametrized using a model based on 
that of D'Ambrosio, Isidori, and Portoles(DIP)^, but with an additional constraint that ensures 
the coupling vanishes at large momenta 1 ^. In terms of the remaining free parameters, the form 
factor is: 

, . 1 -/i(l + a)(xi + x 2 ) 
/dip (xi , x 2 ; a) = — — , (3) 

(1 - (J,Xl)(l - /jLX 2 ) 

where \i = M 2 /M 2 ra 0.032. 

The parity properties of the decay can be extracted from the angle (j) between the planes 
of the two Dalitz pairs in Fig. [2j where pair 1 is defined as having the smaller invariant mass. 
The distribution of this angle from the dominant direct contribution has the form dT/d(j) ~ 
1 — Acos(2(j)) + B sin(2^>), where A « 0.2 cos (2Q) and B m O.2sin(2£)cos<5. A pure pseudoscalar 
coupling, therefore, would produce a negative cos(2(/>) dependence. 

The branching ratio measurement, which we describe here first, makes use of a normalization 
mode in which two pions decay via 7r° — > e + e~7 and the third tt° — » 77. This "double single- 
Dalitz" mode, denoted Kl —> 7r ° 7r D 7r D w here 7r D refers to ir° —* e + e~7, has the same final state 
particles as the signal mode. Both modes are fully reconstructed in the detector and the total 
invariant mass is required to match the kaon's. The two modes are distinguished by a \ 2 formed 
of the three reconstructed 7r° masses. This serves to identify the best pairing of particles for a 
given decay hypothesis, as well as to select the more likely hypothesis of the two. The similarity 
of these modes allows cancellation of most detector-related systematic effects in the branching 
ratio measurement, but also allows each mode to be a background to the other. 

Radiative corrections complicate the definition of the Dalitz decays in general. We define 
the signal mode 7r° — ► e + e~e + e~ to be inclusive of radiative final states where the squared 
ratio of the invariant mass of the four electrons to the neutral pion mass x^ e = (M^/M^o) 2 is 
greater than 0.9, while events with x^ e < 0.9 (approximately 6% of the total rate) are treated 
as 7T° — > e + e~e + e~7. For normalization, the decay 7r° — > e + e _ 7 is understood to include all 
radiative final states, for consistency with previous measurements of this decay Radiative 
corrections in this analysis are taken from an analytic calculation to order 0(a 2 ) ©. 

Radiative corrections complicate the definition of the Dalitz decays in general. The signal 
mode 7T° — > e + e~e + e~ is defined to be inclusive of radiative final states where the squared 
ratio of the invariant mass of the four electrons to the neutral pion mass x^ e = (M^/M^o) 2 is 
greater than 0.9, while events with x^ e < 0.9 (approximately 6% of the total rate) are treated as 
7T° — > e + e~e + e~7. Radiative corrections in this analysis are taken from an analytic calculation 
to order 0(a 2 ) ^. 

The final event sample contains 30 511 signal candidates with 0.6% residual background 
and 141 251 normalization mode candidates with 0.5 % background (determined from the Monte 
Carlo simulation). The background in the signal sample is dominated by mistagged events from 
the normalization mode, v KTeV finds the following the ratio of decay rates: 

B x>0.9 . B 

eeee 21 = 0.2245 ± 0.0014(stat) ± 0.0009(syst) . (4) 

The 7T° — > e + e~e + e~ branching ratio can be calculated from the double ratio using the known 
values £ 77 = 0.9980 ± 0.0003 and B eei = (1.198 ± 0.032) x 10" 2 M This yields B%>™ = 
(3.26 ± 0.18) x 10 -5 , where the error is dominated by the uncertainty in the ir° — > e + e 7 
branching ratio. KTeV uses the radiative corrections mo del© to extrapolate to all radiative final 
states, finding: 

eeee (7) • 77 = 23g3 ± .0015(stat) ± 0.0010(syst) , (5) 
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Figure 3: Distribution of the angle <f>, in units of 7r, between the planes of the two e + e pairs for 7r° — > e + e e + e 
candidate decays. The solid histogram shows the Monte Carlo expectation for negative parity. 



and l? eeee ( 7 ) = (3.46 ± 0.19) x 10 5 . This branching ratio result is in good agreement with 
previous measurements^. 

The parameters of the 7r°7*7* coupling are found by maximizing an unbinned likelihood 
function composed of the differential decay rate in terms of ten phase-space variables. The first 
five are (x\, X2, yi, 1/2, 4>), where x±, X2, and 4> are described above and the remaining variables 
yi and y2 describe the energy asymmetry between the electrons in each Dalitz pair in the ir° 
center of mass^. The remaining five are the same variables, but calculated with the opposite 
choice of e + e~ pairings. The likelihood is calculated from the full matrix element including the 
exchange diagrams and 0(a 2 ) radiative corrections. 

The fit yields the DIP a parameter and the (complex) ratio of the scalar to the pseudoscalar 
coupling. For reasons of fit performance, the parity properties are fit to the equivalent parameters 
k and rj, where K + irj = tan £e 1<5 . The shape of the minimum of the likelihood function indicates 
that the three parameters a, k, and rj are uncorrelated. Acceptance-dependent effects are 
included as a normalization factor calculated from Monte Carlo simulations. 

Systematic error sources on a and k are similar to those for the branching ratio measurement. 
The dominant systematic error is due to variation of cuts, resulting in a total systematic error 
of 0.9 and 0.011 on a and k respectively. For the rj parameter, the primary uncertainty results 
from the resolution on the angle (j) between the two lepton pairs. This behavior was studied with 
Monte Carlo simulation and a correction was calculated. The uncertainty on this correction 
results in a systematic error of 0.031. 

The 4> distribution is shown in Fig. [3l For plotting the data a unique pairing of the four 
electrons is chosen such that x\ < X2 and the product X1X2 is minimized: this choice represents 
the dominant contribution to the matrix element. It is clear that the pseudoscalar coupling 
dominates, as expected, with no evidence for a scalar component. The distributions of all five 
phase space variables agree well with the Monte Carlo simulation. 

The parameters n and 77 are transformed into limits on the pseudoscalar-scalar mixing angle 
£ under two hypotheses. If CPT violation is allowed, then the limit is set by the uncertainties in 
77, resulting in Q < 6.9° at the 90% confidence level. If instead, CPT conservation is enforced, r\ 
must be zero, and the limit derives from the uncertainties on k, resulting in £ < 1.9°, at the same 
confidence level. These limits on £ limit the magnitude of the scalar component of the decay 
amplitude, relative to the pseudoscalar component, to less than 12.1% in the presence of CPT 
violation, and less than 3.3% if CPT is assumed conserved. The limits on scalar contributions 
apply to all 7r° decays with two-photon intermediate or final states. 

This analysis confirms the negative parity of the neutral pion with much higher statistical 
significance than the previous result, and places tight limits on nonstandard scalar and CPT- 
violating contributions to the 7r° — > e + e~e + e~ decay. 



3 Lepton Flavor Violation 



Lepton Flavor Violation (LFV) in weak decays is a key signature of severa l beyond-Standard 
Model physics scenarios. SupersymmetryQJ] new massive gauge bosonsQ^HS an d technicolor^^ 
all can lead to LFV decays which might be within reach of current experiments. Searches in Kl 
decays are complementary to searches in the charged lepton sector, since Kl decays probe the 
s — > dfie transitional. KTeV-E799 has searched for the decays Kj, — > 7r°/i =l= e T and ir° — > fj/^e^, 
and has made the first reported search for Kl — > 

In each case, the analysis required two charged tracks, one of which was identified as a muon 
and the other an electron. The key detector elements for particle identification were E/p in 
the Csl calorimeter, response of the TRD, and muon hodoscopes downstream of the muon filter 
steel. Clusters in the Csl with no tracks pointing to them were considered photons. 



The dominant background for Kl — > 7r /i ± e T was the decay Kl — * vr^e^fe (K e ^), with a -k^ 
decay or punch through to the muon hodoscopes, accompanied by two accidental photons faking 
a 7T°. Since accidental photons were often accompanied by other accidental activity, we removed 
events with evidence of additional in-time activity in the detector. Additionally, the two photons 
were required to form a good 7r° mass, and the square of the 7r° momentum in the Kl rest frame 
was required to be positive and therefore physical. 

The signal and control regions were defined using a likelihood variable L derived from pf , the 
sum of the momentum components of all final-state particles perpendicular to the kaon flight 
line, and M^o^, the invariant mass of the 7r°/xe system. The signal (control) region was defined 
by a cut on L chosen to retain 95% (99%) of signal Monte Carlo events after all other cuts. 
Expected background levels were 0.66 ±0.23 events in the signal region and 4.21 ±0.53 events 
in the control region. Both the signal and control regions were blind during the analysis. Figure 
H] shows the pi — M^o^e plane after all cuts: five events were found in the control region and 
zero in the signal. The resulting limit is B(Kl — ► vr°/i ± e =F ) < 7.56 x 10~ n at 90% CL, a factor 
of 82 improvement over the previous best limit for this mode.^^ 
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Figure 4: Surviving events in the pt — M n o^ e plane for the Kl — > 7r°^i ± e T search data. The signal and control 

regions are shown as the inner and outer solid contours. 



3.2 Other lepton flavor violating modes 



KTeV has also searched for the decay Ki — > Reconstructing a second ir° greatly 

reduces the backgrounds, so some particle identification and anti-accidental cuts were relaxed 
to improve the signal acceptance. A similar analysis, including a cut on a kinematic likelihood 
variable, yielded no events in either the control region or the signal region. This resulted in a 
limit B(Ki — > 7r°7r /i =l= e^) < 1.64 x 10~ 10 . This is the first limit reported for this decay. 

The decay chain Kl — > 7r°7r 7r°, tt° — > [i^e^ gives the same final state particles as Kl — > 
and therefore the same analysis procedure applies with the additional requirement 
that the invariant mass M^ e ~ M^o. Since no events were found, the limit is B(ir° — * /U^e^) < 
3.59 x 10~ 10 . This limit on 7T° — > n is equally sensitive to both charge modes, while the 
previous best limits were not^-l Assuming equal contributions from both charge combinations, 
KTeV's result is about a factor of two better than the previous best limit on ir° — > fj, + e~ and 
about a factor of 10 greater than the previous best limit on ir° — > fi~e + . 

4 Conclusion 

KTeV has completed several measurements recently on the decays of neutral K and ir mesons. 
The measurement of 7r° — > e + e~e + e~ represents the best direct determination of the parity of 
the 7T° and the first searches for nonstandard parity and CPT violation in this mode. It also 
yields the best branching ratio and the first measurement of the form factor in this mode. The 
limits on lepton flavor violation are now the most stringent in the world for these decay modes. 
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